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The kinetics of an irreversible catalytic reaction on a substrate of arbitrary dimension is examined. 
In the limit of infinitesimal reaction rate (reaction-controlled limit), we solve the dimer-dimer surface 
reaction model (or voter model) exactly in arbitrary dimension D. The density of reactive interfaces 
is found to exhibit a power law decay for D < 2 and a slow logarithmic decay in two dimensions. 
We discuss the relevance of these results for the monomer-monomer surface reaction model. 

PACS numbers: 05.40. +j, 68.10.Jy, 82.20.Mj 



In modeling heterogeneous catalysis |jj , the monomer- 
monomer surface reaction model plays an important role, 
at least from the theoretical point of view since an ap- 
pealing simplicity of this model allows one to examine 
several issues analytically. In particular, investigations 
of the monomer- monomer model clarified the role of fluc- 
tuations [p|- |8|, interfacial roughening ||, d iffusion of the 
adsorbants |10|, and surface disorder [|11|. In the sim- 
plest situation (no diffusion, no disorder, etc.), it was 
found that single-species clusters grow with time when 
the dimensionality D of the substrate is sufficiently small, 
D < 2. However, the details of the coarsening like the de- 
cay rate of the density of reactive interfaces remain uncer- 
tain in two dimensions, - simulations |^ |P2|Jilj| ] revealed a 
very slow decay which could be logarithmic or power law 
with a small exponent. In this paper, we clarify these 
questions by computing analytically kinetic character- 
istics of an idealized version of the monomer- monomer 
model, the voter model. We then expand these results 
and perform numerical simulations for the full model. 

The monomer-monomer surface reaction process can 
be schematically represented by the following kinetic 
steps: 
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V^A V , 
V^By, 



(1) 



A v + B V -Z+AB \ + 2V. 



A and B particles impinge upon a surface, with respec- 
tive rates kA and ks, and adsorb onto vacant sites V 
to form a monolayer of adsorbed particles, Ay and By. 
Nearest-neighbor pairs of dissimilar adsorbed particles, 
Ay By, react and desorb with rate k r , leaving behind two 
vacancies. For Ua ^ ks, the adsorption imbalance leads 
to the quick saturation of the surface by the majority 
species. For kA = ks and for dimensions D < 2, there 
is a fluctuation-induced coarsening of the surface into 
growing A and B adsorbed islands. This nontrivial case 
of equal adsorption rates will be considered in the follow- 
ing. Furthermore, in theoretical analysis we will restrict 
ourselves to the reaction- controlled limit, k r <C kA = kg, 



which was found to provide qualitatively the same be- 
havior as the general case [[|J(| but more amenable to 
theoretical treatment. 

In the reaction-controlled limit, the substrate quickly 
becomes completely covered and then stays covered for- 
ever, since in units of the typical, i.e. adsorption, time in- 
terval unoccupied sites are refilled instantaneously. The 
kinetics of the monomer-monomer surface reaction model 
is conveniently described by a mapping onto the Ising 
model with mixed zero-temperature voter dynamics and 
infinite-temperature Kawasaki dynamics ||. Remember 
that in the voter model jl4j , sites have two opinions which 
can be marked by A and B. Each site keeps its opinion 
some time interval, distributed exponentially with char- 
acteristic time t, and then assumes an opinion of a ran- 
domly chosen neighboring site. If a site is surrounding 
by similar sites it does not change its opinion and there- 
fore the voter dynamics is zero-temperature in nature. 
In the original monomer- monomer model, the reacting 
neighboring sites, Ay and By, desorb and then unoccu- 
pied sites immediately refilled by Ay By (no reaction), 
Ay Ay or By By (voter dynamics), or By Ay (Kawasaki 
exchange dynamics whose effective temperature is infi- 
nite since the reaction rate does not depend on the con- 
tent of neighboring sites). Thus the voter model can 
be considered as an idealized variant of the monomer- 
monomer model. 

Interestingly, the voter model can be exactly mapped 
onto the reaction-controlled limit of the dimer-dimer sur- 
face reaction model jp|. Indeed, in the dimer-dimer 
model an empty pair that appears after the reaction 
event, A v +B v -> AB f +2V, is refilled by an A 2 or a B 2 
dimer, so the resulting dynamics is identical to the voter 
model dynamics. However, natural initial states for the 
voter model and the dimer-dimer surface reaction model 
are different. For the voter model, a lattice completely 
covered by A and B monomers without correlations in 
initial positions provides a reasonable initial condition. 
For the dimer-dimer model, if we start from an empty lat- 
tice and fill it by random sequential adsorption of dimers 
one cannot reach a fully covered lattice and instead ap- 
proach a so-called jammed state with some single-cite 
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vacancies |Lq]; this jammed state provides natural initial 
condition for the dimer-dimer model. In the following, 
we will consider only completely covered initial states. 

A remarkable feature of the voter model is its solv- 
ability. It is evident in one dimension where the voter 
model is identical to the kinetic Ising model with zero- 
temperature Glauber dynamics |l4|,[ll|. Surprisingly, the 
voter model can be solved in arbitrary dimension and 
thus appears to be one of a very few models which are 
solvable in any D. Moreover, the voter model can be 
solved on arbitrary lattice although for simplicity we will 
focus on a (hyper)cubic lattice. To reveal the solvabil- 
ity of the voter model, the spin formulation is convenient. 
Identifying A's (B's) with + (— ) spins so that the state of 
the substrate is described by 5 = [5k], k = (fci, . . . , kn), 
one can verify that the spin-flip rate, W k (5) = W(Sk — > 
— 5k), reads 



Wk(S) = i(l-^S k 5> k+e ^ 



(2) 



Here the sum in the right-hand side runs over all 2D near- 
est neighbors and r defines the time scale of the process. 
In the following, we will set r = A/D to simplify numer- 
ical factors in equations for correlators. The probability 
distribution P(5, t) satisfies the master equation 



dt 



P(S, t)=J2 [^ k (5 k )P(5 k , i) - W k (S)P(S, t)] , (3) 



where the state 5 k differs from 5 only at the site k. One 
can then derive a set of differential equations for the spin 
correlation functions 

(5 k ...5i) = ^5 k ...5iP(5,i). (4) 



For the single-body correlation functions we get Q 

4|(5k> = A k (5k>. (5) 
Here Ak denotes a difference Laplace operator, 

A k (5 k ) - -2,0(50 + ^(Sk+eJ- (6) 



For the two-body correlation functions one has 



4 |(5 k 5 1 ) 



(A k + A0(5k5i). 



(7) 



Similar equations can be written for higher-body correla- 
tion functions. An important feature of these equations 
which allows an analytical treatment is their recursive na- 
ture, - to solve for n-body correlation function one does 
not need the higher correlation functions. The struc- 
ture of equations for correlators is similar to the one that 



arises in the one-dimensional kinetic Ising model with 
zero-temperature Glauber dynamics [ fl6| which is equiv- 
alent to the voter model in ID. 

The general solution to Eqs. (||) reads 



<5 k 



_ P ~Dt/2 



£ViI k -i(V2). 



(8) 



Here Ik(x) is the shorthand notation for the multi-index 
Bessel function, Ik(x) = Y\ 1 < J < D h d with /„ being 
the usual modified Bessel function, and a\ — (S\)(t = 0). 

Although the evolution rules of the voter model do 
not preserve locally the densities of A's and B's, Eq. (||) 
shows that 5~J k (5 k ) = Xa CT i' the total densities are con- 
served. Note, however, that for any finite substrate in- 
tuitively more appealing behavior emerges: The effect of 
fluctuations leads to saturation, i. e. all voters eventually 
share the same opinion thereby stopping the dynamics. 

To find the two-body correlation functions, we first 
make a simplifying assumption that the initial state is 
translationally invariant. Then (5 k 5i) will depend only 
on m = k — 1 for t = 0. Clearly, this holds for later 
times, too. In this situation, the shorthand notation 
Rm = (5k5i) will be used. Thus for the translation- 
ally invariant initial conditions, Eqs. (0) simplifies to the 
lattice diffusion equation 



2 Ajn-Rj, 

at 



(9) 



which should be solved subject to the boundary condition 



Ro(t) = l, 



(10) 



since Rq — (5 k ) =1. It is natural to choose an uncorre- 
lated initial state, 



R m (t = 0) = 0, for all m ^ 0. 



(11) 



Eq. (JsJ) is identical to Eq. (||) up to a numerical factor. 
Therefore, if we forget for a moment about the boundary 
condition of Eq. (|i~(i|), we can use Eq. (p|) to get R m = 
e~ Dt I m (t). However, _R = [ e ~*A)(0] disagrees with 
the boundary condition. To remove this discrepancy it 
is useful to consider the initial-value problem, Eqs. (|^)- 
as the problem with a localized constant source at 
the origin and therefore to look for the solution of the 
form 



R m = e 



-Dt 



I m {t)+ I drJ D (t-T)e- DT I m {T). 
Jo 



(12) 



Mathematically, Eq. (|T^) is a linear combination of exact 
solutions to Eq. (|9|) and therefore (|l2|) also solves the 
linear Eq. (^). Physically, Eq. (^) corresponds to the case 
of initial source i?o|t=o = 1 at the origin, supplemented 
by an additional input Jo{r)dT which is added into the 
origin during the time interval (r, t + dr) to keep the 
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overall density at the origin unchanged, Ro(t) = 1. Thus 
the input strength Jo obeys 



l-[e-%(t)] J 



dTj D (t-T)[e- T I (T)] D . (13) 



The Laplace transform of the strength, Jd(A) = 
/ °° dte~ xt Jo{t), can be expressed through the Laplace 
transform To (A) of the function T]j(t) = [e^ t Io(t)] D , 



Jd(A) = -1 + 



ATd(A) 



Making use of the integral representation jlj 



1 

2^ 



c?q e 



i cos q 



(14) 



(15) 



we express To (A) through the so-called Watson integrals, 

/" 27r rf D q 1 
T " (A) = io ^X + D-E^ncosqi (16) 

Combining these finding yields 



in ID and 

J 2 (A) = -1 + 7T 



A + 2 



2A 



A + 2 



(17) 



(18) 



in 2D. In Eq. (Jig), if (a;) is the complete elliptic inte- 
gral of the first kind, ^(x) = f* /2 d6(l - x 2 sin 2 O)- 1 ' 2 . 
The final expressions for Jn(t) are given by the inverse 
Laplace transform. 

From the two-body correlation functions, physically 
interesting quantities can be found. One such quan- 
tity, the concentration CUs(t) of reactive interfaces, or 
nearest-neighbor adsorbed AB pairs, is given by Cab — 
(1 — i?e;)/2. A straightforward computation gives 



Cab (f) 
1 



-Dt T D-1 



+ 2 J dTj D {t-T)e- UT I^-\r)[h{T)~h{T)]. (19) 

The long-time behavior is obtained by analyzing the 
low- A limit. From Eqs. ( |l6| ) and (|lj) we get 

r \-°/ 2 D<2 

Jc(A)~ I X-Hn-^l/X) D = 2 A^O (20) 
I A" 1 D>2 

which imply 



Jd® ~ < (lnt)- 1 D = 2 
I 1 D > 2 



i-s-oo (21) 



Combining Eqs. (Jl£ 
lations 7 (*) — -fi (*) — 

]l7| one finds the asymptotic behavior for the concentra- 
tion of reactive interfaces 



and(21), and the asymptotic re- 
s'/v^, I (t)-h(t) ~ e*/\/8^F 



r i+D/a £, < 2 

(lnt)- 1 L> = 2 
a-bt- D / 2 D>2 



(22) 



Thus in the long-time limit CUb — > when D < 2, i.e. 
the coarsening, takes place for low dimensional substrates 
while for D > 2 single-species domains do not arise. 

In the borderline two dimensional case the coarsening 
occurs but the concentration of reactive interfaces de- 
creases very slowly. In the recent work jn|, the decay 
of Cab(£) in the voter model has been studied numeri- 
cally. Fitting data by power-law and logarithmic forms, 
Cab ~ t~ u and Cab ~ [hi£] _cr , respectively, the effective 
exponents ui sa 0.096 and a w 0.59 have been observed. 
The theoretical asymptotic value CaB(£) computed from 
Eq. ©is 



° Mt) = 21n(i)+ln(256) + ° 



(23) 



and therefore the asymptotic state arises on a very late 
stage which has not been reached in simulations ||,[l2],|l3| . 
In [ p^| , Evans and Ray did simulations for times t < 
1500t. However, trying to fit the exact result ( |2^ ) with 
Cab ~ [lnt] _cr , one can get at the very best the value 
a = 0.72 when t w 1500t. 

Turn now to the monomer- monomer surface reac- 
tion model in the reaction-controlled limit. The 
voter model solution of Eq. (||) is still valid for the 
monomer-monomer model since addition of the infinite- 
temperature Kawasaki dynamics to the voter model dy- 
namics just results in a change of time scale r to r/2 in 
Eq. (§) j§. Eq. (|) for |k — 1| > 1 also maintains its form, 
up to the replacement r by r/2, while for k — 1| = 1, 
Eq. (Q) undergoes more significant change However, 
in the long-time limit the growth of the characteristic 
space scale makes the underline lattice structure less and 
less important and hence the difference with the voter 
model behavior should decrease with time. So one can ex- 
pect that the main difference between the concentration 
CAB(t) of the voter model and the monomer- monomer 
model lies in the time scale difference. This suggests that 
the asymptotic state for the monomer-monomer model is 
established later than for the voter model giving rise to 
a smaller effective exponent a. This hypothesis is com- 
forted by our simulations where, for longer times than in 
previous simulations [p '. 



12 



13y , we have always found an 
effective exponent a closer to unity. 

We have performed numerical simulations for the 
monomer-monomer model on a square lattice of size 
10 3 x 10 3 in time interval t < 10 5 r. Both the size of 
the system and time span of the simulations were signifi- 
cantly bigger than in previous studies |5|,[l^,[l3| . We have 
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found a ps 0.62 which is to be compared with a ps 0.51 for 
times t < 1500r jlj|. We believe that in a truly asymp- 
totic regime the value a = 1 will appear. Moreover, we 
think that our results for the kinetics of the monomer- 
monomer model of catalysis the reaction-controlled limit 
are qualitatively valid for arbitrary reaction rates. In par- 
ticular, for the opposite extreme, i.e. for the adsorption- 
controlled limit (k,A = <C k r ), we have observed 
a ps 0.42 instead of a ps 0.33 § and a ps 0.26 ||. Large 
time simulations are again closer to the true asymptotic 
regime where we think, (hit) -1 behavior should be re- 
covered. 

In summary, for the voter model in arbitrary dimen- 
sion we have found the exact expression for the two-body 
correlation functions. In the most interesting two dimen- 
sional case our exact solution reveals, on the language of 
the catalysis model, that the density of reactive interfaces 
exhibits inverse logarithmic decay. It would be very inter- 
esting to find exact solutions for higher-body correlators. 
Given the fact that one dimensional voter model is com- 
pletely solvable jl8| , exact results for low-body correlators 
in arbitrary dimension can indicate on the complete solv- 
ability of the voter model for any D. Another direction 
of further research is to consider the dimer-dimer model 
of catalysis in the adsorption-controlled limit. Contrary 
to the monomer-monomer model, where there is some 
evidence that the qualitative behavior of the system is 
the same in the reaction- and adsorption-controlled limit, 
the behavior of the dimer-dimer model should be very 
different in these two limiting cases. In the adsorption- 
controlled limit, an infinite number of adsorbing states 
can play a significant role in the dynamics (existence of 
adsorbing states is evident, e.g. in one dimension any 
state with A- and 5-islands separated by single empty 
sites will be an adsorbing state). Note that there are just 
two trivial adsorbing states for the monomer- monomer 
and monomer-dimer models so the dimer-dimer model is 
very different from these previously studied models. Rich 
kinetic behaviors, resembling the ones of the deposition- 
evaporation models |l9|j20|, can be envisioned. 
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